Dislocation core field. I. Modeling in anisotropic linear elasticity theory 
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Aside from the Volterra field, dislocations create a core field, which can be modeled in linear 
anisotropic elasticity theory with force and dislocation dipoles. We derive an expression of the 
elastic energy of a dislocation taking full account of its core field and show that no cross term exists 
between the Volterra and the core fields. We also obtain the contribution of the core field to the 
dislocation interaction energy with an external stress, thus showing that dislocation can interact 
with a pressure. The additional force that derives from this core field contribution is proportional 
to the gradient of the applied stress. Such a supplementary force on dislocations may be important 
in high stress gradient regions, such as close to a crack tip or in a dislocation pile-up. 

PACS numbers: 61.72.Lk, 61.72.Bb 



I. INTRODUCTION 

Far from a dislocation, the elastic field is well described 
by the Volterra solution^ : this predicts that the displace- 
ment and the stress are respectively varying proportion- 
ally to the logarithm and the inverse of the distance to the 
dislocation. But the elastic field may deviate from this 
ideal solution close to the dislocation. Such a deviation 
corresponds to the dislocation core field. It arises from 
anharmonicities in the crystal elastic behavior, especially 
in the high-strained region of the core, as well as from 
perturbations caused by the atomic nature of the core. It 
is, in part, responsible for the dislocation formation vol- 
ume, which manifests itself experimentally^ through an 
increase of the lattice parameter with the dislocation den- 
sity. This leads to an interaction of the dislocation with 
an external pressure^. Although the core field decays 
more rapidly than the Volterra field, it can modify the 
elastic interaction of dislocations with other defects^"^. 
For instance, equilibrium distances in a dislocation pile- 
up are affected by this core field''. As a consequence, 
the stress concentration at the tip of the pile-up is en- 
hanced. This can favor fracture initiation or yielding for 
an edge or mixed dislocation pile-up, or cross-slip for a 
screw pile-up^. The stress produced by this core field also 
tends to open the crack in front of the dislocation pile-up 
on the glide plane®. This explains the nucleation of a 
crack in a mixed mode I-II or I-III sometimes observed 
experimentally in the glide plane of the pile-up. Without 
the core field, only modes II and III would be possible. 
The core field can also alter dislocation properties such 
as their elastic energy^'^ or their dissociation distance in 
fee metals^'^. Finally, it contributes to the elastic inter- 
action between dislocations and impurities, and thus may 
explain part of the solid solution hardening^°. 

One can use atomistic simulations, either based on em- 
pirical potentials or on ab initio calculations, so as to take 
full account of the core field when studying dislocations. 
On the other hand, the core field can be modeled within 
linear elasticity theory, using an equilibrated distribution 
of line forces parallel to the dislocation and located close 
to its core^'^^"^'^. A multipole expansion of the distribu- 



tion leads then to an expression of the core field in term 
of a series. Usually, only the leading term of this series is 
considered. The core field is then equivalent to an ellipti- 
cal line source expansion and is fully characterized by the 
first moments of the line force distribution. Comparison 
with atomistic simulations have shown that this approach 
correctly describes the dislocation core field'*''^'^'^^'^'*. 



Until now, only few studies^'^'^'^ have included this 
core field in the calculation of the elastic energy of dis- 
locations or of their interaction with an external stress 
field. Most of the time, dislocation elastic energies are 
obtained by considering only the Volterra elastic field. 
Such an approximated approach may lead to some errors. 
Notably, simulation boxes used within ab initio calcula- 
tions are usually too small to neglect the core field. A 
recent study of the screw dislocation in iron^ has shown, 
indeed, that it is necessary to include this core field in 
the computation of the elastic energy when deriving core 
energies from ab initio calculations. The purpose of this 
paper is to extend the modeling approach^ ^ of the core 
field within anisotropic linear elasticity so as to include 
its contribution in energy calculations. Previous studies, 
which have considered this core field contribution, either 
assumed that the elastic behavior is isotropic^'^'^'^, or 
that the elastic constants obey a given symmetry^'^^ in- 
compatible for instance with a (111) screw dislocation in 
a cubic crystal. 



In this paper, we first review how the elastic field of a 
line defect, including the core field, is modeled within lin- 
ear anisotropic elasticity theory. The approach of Hirth 
and Lothe^-'^ is generalized so as to describe the core field 
not only through a distribution of line forces but also a 
distribution of dislocations. The elastic energy of the line 
defect is then computed, thus showing the extra contri- 
bution arising from the core field. Finally, we determine 
the infiuence of the core field on the interaction of the 
line defect with an external stress. 
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II. ELASTIC FIELD OF A LINE DEFECT 

We consider a static line defect, the line direction of 
which is denoted 63. Such a defect can be cither a dislo- 
cation, a line force, or the combination of both. Eshelby 
et al}^ have shown that the elastic displacement and the 
stress created at a point of coordinates x can be, respec- 
tively, written as 

1 ® 

Uk{^) = -Y.A'^M, (la) 

a = l 

a,,(x) = l^Bg.,A^^:g^, (lb) 

a=l 

where the variable Za is related to Cartesian coordinates 
by Za = xi + PaX2- The matrix iJ^".^ is obtained from 
elastic constants Cijki by 

^Tjk — C'yfel + PaCijk2- 

The six roots Pa are the imaginary numbers, for which 
the following determinant is null 

|{Bfife +p„i?f2 J|=0. (2) 

The vectors A'^, associated to each root Pa, are the non- 
null vectors that verify the relation 

{Bt,k+Pc.Bt^k)At=0. (3) 

In all the above expressions and in the followings, we use 
the Einstein summation convention on repeated indexes, 
except for Greek indexes a that design the different roots 
Pa- Such summations on a will always be explicitly in- 
dicated as in Eq. (1). 

The six roots pa are necessary complex. If p„ is a 
solution of Eq. (2), then its complex conjugate Pa* is 
also a solution. We sort the roots Pa according to the 
following usual rule 

$>(Pc«) > and Pa+3 = Pa*, Va e [1 : 3], (4) 

where 5(pq) is the imaginary part oi Pa- With such a 
convention, the matrices Bfj,^ verify the relation 

and the vectors A'^ can be chosen so that 

A-+' = Ar. (6) 

As the elastic displacement has to be real, the functions 
fa have also the property 

fa+3{z*) = faizf. (7) 

The general form of the function fa defining the elastic 
displacement and the stress (Eq. 1) is a Laurent series^^. 



If we restrict ourselves to a line defect in an infinite crys- 
tal, the series is limited to the following terms 

■^"(-"^ = ^ 2^ ( + ^ ^'-^ j ' 

\ fc=— oo / 

with i = T. The sign ^ in this equation has to be 
taken as — for 1 < a < 3 (roots having a positive imagi- 
nary part) and + for 4 < a < 6 (roots having a negative 
imaginary part). In (2;) is the principal determination of 
the complex logarithm, whose imaginary part belongs to 
[— TT : 7r[, thus showing a discontinuity in M^. 

Far from the line defect, the main contribution to 
the function fa, and thus to the elastic displacement, 
arises from the logarithm term. This corresponds to the 
Volterra elastic field created by a dislocation and to the 
2D elastic Green function for a line force. Glose to the 
line defect, the other terms in Eq. (8) may also lead 
to a relevant contribution. For a dislocation, these ad- 
ditional terms describe the core field. In the following, 
we only consider the main contribution to the core field 
corresponding to the term fc = 1 in the Laurent series. 
This correctly describes the core field far enough from the 
line defect. The superposition of the Volterra and of the 
core fields, which gives the total elastic field created by a 
line defect, is then obtained from the following truncated 
series 

^"(^^~oo^2^(^«l^(^) + ^«"1)- 

A. Volterra elastic field 

The Volterra elastic field is given by the logarithm in 
Eq. (8). This leads to the following displacement and 
stress fields 

4 W = ^ E T^^S'^^^^";:- (9b) 

Q — 1 

This corresponds to the long range elastic field of a dis- 
location of Burgers vector b or a line force of amplitude 
F if the coefiicients Da verify the system of equations^^ 

1 6 

Y,AtDa = -hk, 

\j2B?2kAlDo. = -Fi. (10) 

Stroh^^'^'' proposed a simple solution to this system of 
equations. In that purpose, he defined a new vector 

Lt = Bt2kAt = -^Bf^^Al. (11) 

Pa 
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As the vectors A" arc defined through the equation (3), 
their norm is not fixed. One can therefore choose their 
norm so that 



2AfLf = l, Va. 



(12) 

Stroh showed that such a definition of the vectors Af and 
Lf leads to the orthogonahty property 



where we have used Eq. (14) to ehminate the first term 
of the expansion. Using Eq. (9) and taking the limit 
a"^ 0, one finally obtains the expression of the core 
field 



1^1 1 
(x) = y Ttt^A'^C-^ 



2^ 2ni 

a=l 



Xi + PaX2 



(15a) 



AfL^ + A^L 



where is the Kronecker symbol. These vectors also 
verify the following relations''^'''^^ 



2 ^ 27ri 

a=l 



with 



J2AfA'}=0, $:^?^i=0 5^A?L," = 5,,. C-' = 2AnMa+PaM,2) + 2LnP^^ 



PaPi2) 



a=l 



These orthogonality properties lead to the expression of 
the coefficient 



D^ = -2mbi+AfFi). 



(13) 



Core field 



The Volterra solution models the elastic field created 
by a dislocation far enough from the dislocation core. 
Close to the core, the dislocation core field may be rele- 
vant too. Wc model this additional clastic field by con- 
sidering the term 1/z in Eq. (8). Gehlen et al.^ have 
shown that this term may be obtained from dipoles of 
line forces. It is also possible to consider dipoles of dis- 
locations, which may be more natural to model the core 
field of a dissociated dislocation. Therefore, we assume 
that the core field can be modeled by an equilibrated dis- 
tribution of dislocations and line forces of force amplitude 
F'^ and of Burgers vector located at a'' . All line force 
and dislocation directions are assumed to be coUinear to 
63. As the distribution is equilibrated, the resultant of 
the forces and the total Burgers vector have to vanish 



^F« = and ^b« = 0. 



(14) 



The elastic displacement and the stress of this distribu- 
tion is given by the superposition of the Volterra elastic 
field created by each line defect 

n^,(x) = 5^«r'(x-a«), 

9 



4.(x) = ^a,f' (x-a"^). 



We then assume that the norm of x is large comparcxl to 
the norm of the vectors a'. One can thus make a limited 
expansion^ 



.9,11,18,19 leading to 



4(x)=-E^C^«^+o(iia'n, 



aa^'"' fx) 



where M and P are respectively the first moment tensors 
of the line force and of the dislocation distribution 



Mi, 



Y^Ffa] and Pi, =Y.^hV 



As we assume that the distribution of line defects rep- 
resentative of the core field is equilibrated, it does not 
produce any torque. This implies that the tensor Mij 
is symmetric^*. The first moment tensors M and P can 
be simply deduced from the homogeneous stress com- 
puted in atomistic simulations using periodic boundary 
conditions^. Another method based on path-independent 
interaction integrals computed through the field observed 
in atomistic simulations has also been proposed^"'^^. 



III. ELASTIC ENERGY OF AN ISOLATED 
LINE DEFECT 

We now calculate the elastic energy of a line defect, 
such as a dislocation, taking into account its core field. 
The elastic field created by the line defect is thus the su- 
perposition of the Volterra solution given by Eq. (9) and 
of the core field given by Eq. (15). We define the elas- 
tic energy of the line defect as the integral of the elastic 
energy density over the volume in-between two cylinders 
centered on the line defect. The inner cylinder of radius 
fc isolates the line defect core: elastic fields are diverging 
at the line defect position and one needs to exclude the 
core region, where elasticity breaks down. The external 
cylinder of radius i?oo is introduced to prevent the elastic 
energy from diverging. Then, Gauss theorem allows us 
to obtain the elastic energy 

E=\l {al+a<r.){v:( + ut)AS,, (16) 

" s 

where the integration surface S is composed of both cylin- 
der surfaces and the branch cut, which isolates the dis- 
placement discontinuity. We consider cylinders of unit 
height so as to express the elastic energy per unit length 
of line defect. 

This elastic energy can be decomposed into three dif- 
ferent contributions: the contribution of the Volterra so- 
lution, the contribution of the core field and the cross 
interaction between both elastic fields. 
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Volterra contribution 



Volterra - core field interaction 



The Volterra contribution corresponds to the prod- 
uct o-JjU^ in Eq. (16). It is given by the well-known 

resulti6.i7,22 

= I {b,K^,bj + F,K':^Fj) In , (17) 

where we have defined the second rank tensors 

6 6 

K-i = Y ^Tr^^^ and K'^. = V Ttt^A'^A'^. (18) 



a=l 



Then we calculate the interaction energy between both 

clastic fields created by the line defect. Two different 
integrals can be used to obtain this interaction energy^*: 



We rather use the first definition to evaluate E^~^: as 
the core field displacement u'^ does not show any discon- 
tinuity except at the origin, the integration surface of the 
first integral is simply composed of the inner and external 
cylinders. This leads to the following interaction energy 



B. Core field contribution 

The contribution of the core field to the elastic energy 
corresponds to the product cr'^jUl in Eq. (16). As the core 
field does not create any displacement discontinuity, the 
integration surface is simply composed of the inner and 
the external cylinders. This leads to the contribution 



1 ^ 1 
8 27ri ' " 



a=l 
1 



The integrals /^(p, q) and Iy{p. q) are defined in the Ap- 
pendix (Eq. (A.l)). We use the fact that I^{p,q) = 
—ply{p,q) as well as the property verified by the vec- 
tors (Eq. 3) and the definition of the vectors Lf (Eq. 
11) to obtain 



a=l 



Iy{Pa,P0) 



Finally, the expression of the integral ly {p, q) given in the 
appendix allows us to write 



E'= 



/ 3 3 



1+Po 



^1 3=1 '^P" 



3*) 



-1* 



■ (19) 



The dependence of this expression with shows that 
the elastic energy of the core field is concentrated close to 
the line defect. It is possible to take the limit i?oo oo, 
and thus to define an elastic energy associated with the 
core field in the whole volume excluding the core region, 
where elasticity breaks down. 



a=l 

6 

E {BLl^iPa,P0) + B^^Jl{p^,pp)) AlDf, 

1 1 

The integrals /J (p, q) and ly {p, q) are defined in the ap- 
pendix (Eq. (A. 2)). As these integrals vanish for any 
values of p and q, this leads to E^^'^ = 0. As a re- 
sult, there is no interaction energy between the Volterra 
elastic field and the core field of the line defect, and the 
elastic energy of a line defect is simply the sum of the 
elastic energies of the Volterra field and of the core field. 
Of course, this is true only when the Volterra and the 
core fields are centered at the same point. This may be 
imposed by symmetry, as for the (111) screw dislocation 
in a cubic crystal. ^'^^ When the Volterra and the core 
fields have different centers, ^■'''^'^^ an interaction energy 
between both elastic fields exists. Such a cross term can 
be simply calculated by considering the interaction of the 
core field with the stress created by the Volterra field, as 
described in the next section. 



IV. INTERACTION WITH A STRESS FIELD 

We now consider the interaction energy between an 
external stress field aff^ and a line defect. The external 
stress can be an applied stress or the stress originating 
from another defect. The line defect is located at the ori- 
gin and its line direction is e^. It is characterized by its 
Burgers vector b, its force resultant F and the first mo- 
ments tensors and Pij. The interaction energy can 
be decomposed into two contributions: the interaction 
with the Volterra elastic field and the interaction with 
the core field. The first contribution is well known. ^'^^ 
For a dislocation, it is given by the integral of a°j^^bidSj 
along the dislocation cut, where dSj is an infinitesimal 
surface vector. For a line force, it is given by the scalar 
product Fiuf^'^, where uf^'^ is the displacement associ- 
ated to the external stress field. We now determine the 
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contribution of the core field to the interaction energy 



A. Core field contribution 

The interaction energy of the core field with the stress 
field crfj^* can be obtained by considering the line defect 
distribution responsible for the core field, thus using the 
same approach as used by Siems^^'^"^ for a point defect. 
The interaction energy is then given by 

q Jo 

where Cjki is the permutation tensor. The first term rep- 
resents the interaction with the different dislocations of 
the distribution (e^fesa^dA is the infinitesimal surface vec- 
tor along the dislocation cut), and the second term rep- 
resents the interaction with the line forces. A limited 
expansion of cr|j'* and of uf^^ at the origin leads to 

q ^ 

We finally use the fact that the tensor Mij is symmet- 
ric and take the limit a' ^ to obtain the interaction 
energy 

winter ^ ^ext^Q) ^^.^^p.^ _ ^..^^Mfc,) , (20) 

where the clastic compliances Sijki are the inverse of the 
elastic constants {SijkiCkimn = 5 {^im^jn + 

Thus Eq. (20) shows that an additional contribution 
needs to be considered in the interaction energy of a 
line defect with a stress when a core field is present. 
In particular, this contribution of the core field leads to 
a dislocation-pressure interaction which can modify the 
kink formation energy and the dislocation line tension at 
high pressures^. Such a dependence of the dislocation 
energy with the pressure has been observed in atomistic 
simulations^^"^^. Eq. (20) should allow to model this de- 
pendence, or at least the first order variation. 

B. Force acting on a line defect 

The external stress field aff * creates on the line defect a 
force which derives from the interaction energy. Without 



the core field, this force would be simply given by the 
Peach-Koehler formula for a pure dislocation and by the 
product —Sijki<Tl^^Fj for a pure line force. Because of 
the core field, there is an additional force acting on 
the line defect. This force derives from the interaction 
energy given by Eq. (20) 

da^f(0) 

fn = {ijksPik — SijklMkl) ■ (21) 

Because of the core field, there is a force acting on the 
line defect, which is proportional to the gradient of the 
applied stress. 



C. Elastic energy of an isolated dipole 

The calculation of the elastic energy of an isolated dis- 
location dipole is one important application, where one 
needs to calculate the interaction energy of a line de- 
fect with a stress field. In that case, the external stress 
field is created by the other line defect composing the 
dipole. Here we determine the elastic energy of a dis- 
location dipole that is assumed to be isolated from any 
other defect. This elastic energy is defined as the inte- 
gral of the energy density on the whole volume except 
two cylinders of radius Vc excluding the regions around 
the dislocation core. As the elastic energy is now converg- 
ing, we do not need to introduce an external cylinder as 
we did for an isolated dislocation. The dipole is com- 
posed of two line defects of opposite Burgers vector b 
and opposite force resultant F having the same core field 
characterized by the moment tensors Mij and Pij. We 
assume that 63 corresponds to the line direction of the 
dislocations. The dipole is then defined by the distance d 
between the two dislocations and by the angle (j) between 
the dipole direction and a reference vector ei. 

If the elastic field created by each dislocation compos- 
ing the dipole is only of the Volterra type [Eq. (9)], the 
elastic energy of the dipole is 

El^oie = {biK°jbj + F.K'^.Fj) In (^^) + 2i?7 (</.), (22) 

where the tensors Kfj and K'^j are given by Eq. (18). 

is the contribution from core tractions to the elastic 
energy. Such a contribution arises from the work done 
by the tractions of the Volterra elastic field exerted on 
the cylinders that isolate the dislocation cores. It exists 
even when the core field is neglected and it is given by^^ 
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3 6 



Q=l 8=4 



1 1 

+ 2 E ± 2^ i^i^'mh - Pi^t^^Pj) In (cos <^ + Pa Sin 0). 

c«=l 



Considering now that a core field as described by 
Eq. (15) is also created by each dislocation, one has to 
add to the clastic energy of the dipolc [Eq. (22)] the con- 
tribution from the core field of each dislocation, 2E'^, as 
given by Eq. (19) in the limit i?oo oo. 

Another contribution also needs to be taken into ac- 
count in the elastic energy when dislocations create both 
a Volterra and a core field. It arises from the interaction 
of the total stress field created by the first dislocation 
with the core field of the second one, and vice versa. 
This interaction energy can be calculated using Eq. (20), 
which leads to the result 

^^dipo'ie = (2^.^(d) + 4-(d)) {ejksPik - SijklMki) , (23) 

where the vector d is defined by the coordinates 
d(coS(/), sin0, 0). Equation (23) shows that the elastic 
energy of the dipole now contains a contribution varying 
with the inverse of the distance d and another contribu- 
tion varying with the square of the inverse of d. 

D. Dislocation dipole in periodic boundary 
conditions 

When stiidying dislocations in atomistic simulations, 
one can use periodic boundary conditions. A dipole is 
introduced to ensure that the total Burgers vector of the 
simulation box is null. Atomic simulations allow obtain- 
ing the excess energy associated with the defects present 
in the simulation box. One can deduce from this quan- 
tity dislocation intrinsic energy properties such as their 
core energy. To do so, one needs to calculate the elas- 
tic energy contained in the simulation box. This elastic 
energy includes the elastic energy of the primary dipole 
as well as half the interaction energy with all its peri- 
odic images. When the simulation box is small, as in 
ab initio calculations, one needs to consider not only the 
Volterra field but also the core field of the dislocations 
when computing the elastic energy^. 

The elastic energy of the primary dipole is given in 
the preceding section. The interaction energy between 
two dipoles can be obtained by decomposing it into the 
contributions arising from the different constituents of 
the elastic field. The interaction energy arising from the 
Volterra field of each dipole is obtained thanks to the 
expression given by Stroh^^ for the interaction energy 
between two dislocations. If the coordinates of the vec- 



tors joining the two dislocations are {xi,X2,X3), this part 
of the interaction energy is given by 

^int-eT = - E ^^b^^^t^pf In [X^ +p^X2), 
a=l 

where b^^^ and b^^^ are the respective Burgers vectors of 
each dislocation^''. 

The part of the interaction energy arising from the 
core field is obtained thanks to Eq. (20). The external 
stress crfj^* appearing in this equation corresponds to the 
stress created by the other dislocations, where both the 
Volterra and the core fields are considered. 

When summing all contributions from the different pe- 
riodic images, one should be aware that the sums arc only 
conditionally convergent. This convergence problem can 
be easily resolved using the regularization method of Cai 
et aL^°. 



V. CONCLUSIONS 

We have extended in this paper the approach of Hirth 
and Lothc^^ to model dislocation core fields within lin- 
ear anisotropic elasticity theory by deriving the elastic 
energy of a straight dislocation while taking full account 
of its core field. The obtained expression shows that this 
elastic energy is the sum of the energies corresponding 
to the Volterra field and to the core field, and that no 
cross interaction exists between these two elastic fields. 
We have also shown that the core field leads to an addi- 
tional contribution to the interaction energy between a 
dislocation and an external stress. Through this contri- 
bution, the energy of a dislocation can depend on the ap- 
plied pressure. This interaction with the core field is also 
responsible for an additional force acting on the disloca- 
tion, which is proportional to the gradient of the applied 
stress. Dislocation properties may therefore be affected 
in regions where a high-stress gradient is present such as 
in a dislocation pile-up'''^ or close to a crack tip. 

The interaction of the dislocations caused by their core 
field is shorter-range than their interaction through their 
Volterra field. It will therefore affect the interaction be- 
tween dislocations when they get close enough. Such a 
situation may arise in atomistic simulations, where the 
size of the simulation box may be too small to neglect the 
influence of the core field. ^ One should then take account 
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of the dislocation core field when calculating elastic ener- to the result 
gies, which can be done using the different expressions of 

this paper. An example is given in the following paper, ^'^ /^(p, q) = Q if > and > 0, 

where the developed formalism is applied to the (111) 

screw dislocation in a-iron. 47rip 
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Appendix: Integrals 



The elastic energy of the core field (c/. §111 B) makes 
the two following integrals appear 



J —n 



(cos^ -|-psin^)(cos^ + g'sin^)^ ' 



sm ( 



(cos 9 +psiii 0) (cos + g sin 0)2 



de. 

(A.1) 



These integrals of a rational function of cos [9) and sin [9] 
can be evaluated using the residues theorem'^^ . This leads 



-47ri 

{p - qf 

47ri 

{p - 





if 3(p) > and 3(g) > 0, 
if 3(p) > and 3(g) < 0, 
if 3(p) < and 3(g) > 0, 
if 3(p) < and 3(g) < 0. 



The two integrals appearing in the interaction energy 
between the Volterra and the core fields of a line defect 
(c/. §111 C) are 



il{p^ q) 



cos 6 



72 

V 



{p,q)= r 

J —TV 



^ (cos 6 +psm. 9) (cos 9 + qsm 6) 
sin 9 



^0 (p^ 2) 

(cos 6' + p sin 6') (cos + g sin 6') 
The residues theorem leads to the result I^{p,q) = 

l^{p,q) = 0- 
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